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We present results from a comprehensive partial-wave analysis of 7r ± p elastic scattering and 
charge-exchange data, covering the region from threshold to 2.1 GeV in the lab pion kinetic energy, 
employing a coupled-channel formalism to simultaneously fit ir~p — > r\n data to 0.8 GeV. Our main 
result, solution FA02, utilizes a complete set of forward and fixed-t dispersion relation constraints, 
from threshold to 1 GeV, and from t = to —0.4 (GeV/c) 2 , applied to the 7rN elastic amplitude. 
A large number of systematic checks have been performed, including fits with no charge-exchange 
data and other database changes, fits with few or no dispersion relation constraints, and changes to 
the Coulomb correction scheme. We have also reexamined methods used to extract Breit-Wigner 
resonance parameters. The quality of fit to both data and dispersion relation constraints is superior 
to our earlier work. The results of these analyses are compared with previous solutions in terms of 
their resonance spectra and preferred values for couplings and low-energy parameters, including the 
nNN coupling constant. 



PACS numbers: 14.20.Gk, 13.30.Ek, 13.75.Gx, 11.80.Et 



. I. INTRODUCTION 

The determination of resonance properties for all accessible baryon states is a central objective in Nuclear Physics. 
Pole positions, decay channels, and branching ratios have been extracted from hadronic reaction data. Further results 
•'■h . for magnetic moments and photo-decay amplitudes have been obtained from real and virtual photon interactions. 
This body of information has been used to test QCD-inspired models and, more recently, lattice calculations. 

Most N and A resonances, listed as 3- and 4-star states in the Review of Particle Properties (RPP) 0, have 
had their existence, masses, and widths determined through single-channel fits to scattering data, with 7rN elastic 
scattering being the predominant source. Similarly, most photo-decay amplitudes have been determined in fits to 
single-pion photoproduction. However, a number of recent studies of eta photoproduction and electroproduction 
data have claimed resonance properties for the N(1535) and N(1520) at variance with these long-standing pion 
production results. Multi-channel analyses have provided further estimates [j| which also tend to disagree with the 
single-channel values. 

These problems have motivated an improved analysis of 7rN elastic scattering in the rjN threshold region. Here, 
we report the results of a coupled-channel fit to 7rN elastic and 77N production data, the ttN — > rjN data having 
been incorporated through a Chew-Mandelstam K-matrix formalism Compared to our previous work (solution 
SM95 |5|), the present FA02 analysis has fitted a larger database, in particular new high-quality data in the A(1232) 
resonance region (described in Section [Hj. In addition to the forward C ± (w) and fixed-t B ± (i/, t) dispersion relation 
constraints already used in SM95, our new solution has been constrained by further forward derivative E (u;) and 
fixed-t C =t (t', t) dispersion relations over broader ranges of energy and four-momentum transfer. 

In Section IlIII we will explain how the dispersion relation constraints were implemented. For comparison purposes, 
fits to the available data were also performed with fewer, and no, dispersion relation constraints. This has helped to 
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gauge the relative effect of such constraints on the quality of fit to data. Many other fits were performed to study 
systematic effects on resonance and dispersion relation parameters, in particular the ttNN coupling constant. We have 
performed fits excluding charge- exchange data, employing different datasets in the A(1232) region, using different 
Coulomb correction schemes, and retaining the older Karlsruhe 7] D-waves at low energies, to gauge the influence of 
these poorly known amplitudes. 

In Section IIV Al we compare the present FA02 solution with both our previous SM95 solution and results from 
previous analyses, as well as giving values for the partial-wave amplitudes. The baryon spectrum, and associated 
couplings are given in Section IIVBI whereas results for the 7rNN coupling constant and other dispersion relation 
parameters are discussed in Section IIV Dl In Section we provide a brief summary and consider what extensions 
of this work can be expected in the future. 



II. DATABASE 

The evolution of our database is summarized in Table [fl Over the course of four previous pion-nucleon analyses 0, 
IE IE Hi our energy range has been extended from 1.1 to 2.1 GeV. Additions to the current database 01 are due mainly 
to measurements in the low (below 300 MeV) energy region. We have also incorporated 173 r\ production data (from 
the process 7r~p — > r/n). Below, we list the recent (post-1995) additions for elastic and charge-exchange scattering. 
The eta-production database, being recently added, contains both new and old measurements. Here, we should note 
that the full database contains conflicting results. Some of these with very large x 2 have been excluded from our fits. 
However, all available data have been retained (the excluded data labeled as "flagged" so that comparisons can 
be made through our on-line facility [Tol ]. Some individual data points were also removed from the analysis in order 
to resolve conflicts or upon authors' requests (since our previous analysis [j|, we have flagged 48 n + and 35 tt~ data, 
most of these being total cross sections). Some of the data, listed as new, were available in unpublished form at the 
time of our previous analys is [5j . A complete description of the database and those data not included in our fits is 
available from the authors [l(j. The energy-angle distribution of recent (post-1995) n + p, n~p elastic, 7r~p — * 7r°n, 
and all ir~p — > r\n data is given in Fig. ^ 

Since most of the new data jH QHl [llllllllllllMIIlllllMIMIlllllllllllare from high-intensity 
facilities (such as LAMPF, TRIUMF, and PSI), they generally have smaller statistical errors and, thus, have a greater 
influence on the fits. As mentioned above, a large fraction of the more recent 7r ± p data were produced at energies 
spanning the A(1232) resonance. These data have been taken mainly at TRIUMF. From this source, we have added 
106 7r + p and 54 n~p differential cross sections, from 140 to 270 MeV at medium scattering angles [Lj, and 41 n + p 
(37 tt~p) cross sections, from 90 to 140 MeV 14] . A further 194 n + p and 81 -K~p A y data between 90 and 280 MeV, 
at medium and backward scattering angles, have been collected using the CHAOS facility at TRIUMF [3 [3. At 
low energies, 55 to 140 MeV, and backward scattering angles, 13 ir + p and 89 n~p A y data have also become available 
from CHAOS [13 ■ 

A few LAMPF experiments have now been analyzed and added to our database. These include 36 high-quality 
polarized charge-exchange data between 100 and 210 MeV 0], 40 low-energy (10 to 40 MeV) cross sections [20|. and 
6 forward charge- exchange differential cross sections at 27 MeV p^ . 

We have added 23 ir + p and 6 n^p TRIUMF [l]| and 44 n+p and 15 n^p LAMPF ^ partial total cross sections 
from 40 to 280 MeV 0], though this quantity is not fitted in the analysis. 

PSI experiments have provided 5 n~p differential cross sections in the backward hemisphere between 45 and 
70 MeV [2i| and 11 A y data at 160 and 240 MeV [24[. After a revised analysis, the Karlsruhe group has pub- 
lished a final set of both 7r ± p low-energy differential cross sections [2^ and analyzing powers |2(|. 

A set of polarization parameters P, R, and A for both n + p and ir~p from 870 to 1490 MeV were contributed by 
the ITEP-PNPI Collaboration [HHg. 

In the A(1232) resonance region, there are two sets of data which are in disagreement, primarily on the rising 
side of the resonance: the first includes the PSI total cross section data of Pedroni et al. |29j and the differential 
cross sections of Brack et al. [l4l [3(j > while a second set includes the total cross sections of Carter et al. and 
the differential cross sections of Bussey et al. [12 . Since the resonant A(1232) (P33) amplitude dominates most 
dispersion relations, discrepant data bases in this energy region are no small concern. In our previous solution p|, 
the Bussey et al. data were assigned 5% normalization uncertainties, and the Carter et al. data were assigned 1.5% 
uncertainties, to resolve a discrepancy between the energy-dependent and single-energy solutions. Subsequent to that 
solution, the normalization uncertainties in the Bussey et al. and Carter et al. data were reevaluated (to 1%), and 
the lowest energy measurements in each of the 7r + p and n~p total cross sections were removed, by one of the principal 
authors [23 . The data were so employed in the present solution. However, with the addition of new TRIUMF 
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differential cross section |l3j and polarization data |l5l Il6j , the changes in phases and dispersion relation parameters 
arising from selecting one of the aforementioned discrepant data sets have greatly diminished, as will be discussed in 
Sec. II V Dl All of the above data sets are included in the fitted database. 

Most measurements of the ir~p — > r\n reaction cross section are rather old and sometimes conflicting. There are 
few cross section measurements above 800 MeV and no polarized measurements below 1040 McV We have 

selected 138 differential and 35 total cross sections for our anal ysis , with an additional 176 data points added to 
the on-line database but not included in this study [M HI MlH HI IME1 El El El El El El A detailed 
analysis of the older data can be found in the review by Clajus and Nefkens |47| . Previous unpolarized measurements 
are listed in Ref. 48] . One new set [34| of data has been included in our anal ysis . Total and differential cross 
sections for 7r~p — * rjn, near threshold, were measured using the BNL-AGS Facility [3J|. Cross sections were obtained 
from threshold (560 MeV) to ~650 MeV with ^10% statistics for the total cross sections. The angular coverage 
was 45—155°. Systematic uncertainties were estimated at ^3%, a value significantly less than claimed in previous 
measurements, where systematics were above 10%. 

III. FORMALISM 

There are three principal components to the methodology we use in fitting nN elastic scattering, charge exchange 
and Tr~p — > r\n data: (i) the energy-dependent parametrization of partial waves, (ii) dispersion-relation constraints, 
and (iii) Coulomb interaction effects. 

A. Chew-Mandelstam K-matrix 

Our energy-dependent partial-wave fits are parameterized in terms of a coupled-channel Chew-Mandelstam K- 
matrix, as described in Refs. 0, . This choice determines the way we introduce and modify the energy-dependence 
and account for unitarity in our fits. 

The scattering into different channels is represented by a matrix T, parameterized in terms of K 1 a (4x4) real 
symmetric K-matrix for each partial wave 

f = K(l-CK)- 1 , (1) 

C being a diagonal Chew-Mandelstam function, with ImCi = pi giving the phase-space function for the i th channel 
(ttN, ryN, irA, and pN). The T-matrix for elastic irN scattering is then given by T v „ = T^^p^; for ttN — > r/N, the 

_ 1/2 

relation being T nr) = ^p^p^) 1 . The pN channel is new to the present solution. Its inclusion has reduced the 
number of parameters by 6, compared to SM95 H, while improving the fits and dispersion relation consistency (the 
p width is accounted for in this approach). It should be noted that the above it A and pN channels are generic and 
included to preserve unitarity. Unlike the SM95 fit, however, the r/N channel has been fitted to physical cross sections. 

In order to control the behavior of each T-matrix near threshold, the K-matrix elements were expanded as polyno- 
mials in the energy variable Z = (W — Wth), where W and "Wth are the center-of-mass (y/s) and threshold energies, 
respectively. Multiplying by an added factor of Z allowed the fixing of scattering lengths through the value of the 
leading term. 

Single-energy analyses were parametrized as 

S e = (l + 2iT e ) =cos(p) e 2iS , (2) 

with the phase parameters S and p expanded as linear functions around the analysis energy, and with a slope (energy 
derivative) fixed by the energy-dependent solution. It should be emphasized that the single-energy solutions are 
generated mainly to search for missing structure. The scatter around an energy-dependent solution is also useful as 
a measure of the uncertainty in a partial wave. The variation from point-to-point should not be taken as proof of 
sharp structures, as these binned fits are far less constrained than the global fits (they do not, for example, satisfy 
the dispersion relation constraints) El]. 

Threshold behavior was determined in the following manner. The S-wave scattering lengths were linked to our 
dispersion relation constraints, as described below. The P-wave scattering volumes were searched without constraint, 
but it was ensured that near threshold these waves followed the appropriate Chew-Low forms [HI which are 
approximations to the respective partial-wave dispersion relations. The D-waves were weakly constrained to the older 
Karlsruhe analysis by assigning 5% errors to the Koch values jjj, and the low-energy behavior of higher waves was 
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fixed to Koch's result Once an appropriate hadronic amplitude was determined, charge corrections were applied 
as described in Sec. IIII CI 

B. Dispersion Relation Constraints 

In general, a fit of the K-matrix elements, expanded in terms of an energy variable, may not result in a form 
satisfying all of the requirements imposed by analyticity and crossing symmetry In our analysis, those requirements 
are addressed at fixed four-momentum transfer t by a complete set of fixed-t dispersion relations (DRs), which are 
handled iteratively with the data fitting, as has been described in Ref. @. The DRs contain subtraction constants 
which should be independent of energy (but which can be functions of four- momentum transfer). After each data- 
fitting iteration, these constants are calculated as a function of energy. \ 2 deviations from the average, at a series of 
energies, are then calculated and included as pseudo-data. Figure |21 shows an example of the subtraction constants 
in the forward C ± (c i j) DRs up to 1 GeV, uj being the pion lab energy. The partial- wave amplitudes and the real 
parts of the DR invariant amplitudes are then adjusted to minimize the x 2 from the sum of data and pseudo-data. 
Compatibility with the DR constraints can be controlled through the errors assigned to pseudo-data in the fit. 

In our previous published analysis we employed constraints from fixed-t DRs for the B±(^, t), v equal to 
(s — u)/4M, and forward C ± (w), invariant amplitudes |53 |. Constraints on the partial- wave fits were generated at 
fixed values of t ranging from to —0.3 (GeV/c) 2 and pion lab kinetic energies from 25 to 600 MeV. The B±(v, t) 
DRs were constructed in the Hiiper form |6j. The full 7rN amplitude was parameterized bytwo invariant amplitudes 
B ± (^, t) and C ± (^, t) (or equivalently B ± (j/, t) and A ± (z/, t)). While in the analysis of Ref. [f| the C* amplitudes were 
reasonably well behaved, improvement called for their inclusion in DR constraints at non-forward angles. 

With the inclusion of constraints from the fixed-t G ± (v, t) DRs, we complete the fixcd-t DRs coverage of the 
full 7rN amp litude. The C ± DRs were chosen since parameterization^ of the very high energy amplitudes are readily 
available [6|. The energy and four-momentum transfer coverage was increased for most fixed-t DRs to < 1000 MeV 
and < — t < 0.4 (GeV/c) 2 . This t range covers the entire angular range up to ~450 MeV, across the important 
A(1232) resonance region, which contains the most accurate and complete data sets. Due to use of analytically 
smooth energy-dependent amplitudes, the DRs are approximately satisfied over wider constraint energy and t ranges 
(see results in Section riVD|) . 

In addition to the fixed-t C ± (^, t) DRs, the present analysis includes constraints from the forward derivative DRs 

E ± M = | [A± (u, t) + tuB± {u, t)]\t=o (3) 

utilized in Ref. [53J- These DRs were also applied from 10 < IV < 1000 MeV. The E ± (u) DRs provide a stronger 
constraint on higher partial waves, which are more poorly known at low energies. Nevertheless, the I — 1 A(1232) 
resonance gives a dominant contribution, as is the case for most of the other DRs we employ. Consistency of the fixed-t 
C (v, t) DRs is illustrated in Fig. [3J As in Fig. |3 the plotted constant is obtained through a delicate cancellation 
between pole contributions and principal value integrals. The variation of these contributions is displayed for the 
forward E ± (w) DRs in Fig. 

Dispersion relations contain principal value integrals over the imaginary parts of invariant amplitudes up to infinite 
energy. The subtraction constants are used to ensure convergence of these integrals. As our analysis extends up to 
2100 MeV, the amplitudes above this energy are obtained from other sources. For the energy region 2100—5000 MeV, 
the Karlsruhe KA84 7rN analysis solution 54] is utilized. At higher energies, parameterizations were taken from 
Ref. 0. These high-energy contributions are relatively small for the Hiiper DR and isovector amplitudes, but are 
important for the isoscalar E + and C + amplitudes. In these latter cases, while high-energy contributions were 
necessary to achieve sufficient convergence, the various forms of parameterization |(| yielded negligible differences in 
the DR constraints. 

As mentioned above, DRs contain a number of parameters which are a priori unknown, namely the various sub- 
traction constants (if inserted for convergence of the DR integrals) and the irNN coupling constant g 2 /An (when the 
Born term appears). These constants are considered free parameters to be determined from the fits to the data and 
the DR constraints. In principle, these subtraction constants can be determined at any energy. We fix them at the 
physical threshold for the forward DRs, and at the unphysical points C ± (0,t) for the fixed-t DRs. The former 
can be expressed in terms of the S-wave scattering lengths (C ± DRs) and the P-wave scattering volumes af + (E ± 
DRs). 

These parameters could be determined in a number of ways. In the present study, we settled on the following 
procedure: for a particular DR and data minimization run, the Tr~p S-wave scattering length a_ was fixed around 
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precise results extracted from the PSI pionic-atom experiment |55|. and the coupling constant g 2 /4ir and scattering 
volume cii + were fixed at chosen values. In essence, the n~p S-wave scattering length a_ is used as an additional 
data point at threshold. The other parameters (including the isovector S-wave scattering length a~) were allowed to 
float freely. As a result, the subtraction constant was fixed for the forward C + and E + DRs, with real parts of the 
amplitudes forced to particular values in the analysis, whereas analogous constants for the other DRs were allowed to 
find their own level. Once a solution with minimal \ 2 was determined, further runs were performed with new sets of 
fixed (i? 2 /47r, af + ) values, producing a mapping over a grid of parameters. The solution FA02 was selected from this 
grid of solutions based on overall x 2 ■ 

Generally, the \ 2 variation of parameters is parabolic in the region near a minimum. The advantage of a grid method 
is that the depth of the (g 2 /Air, af + ) minimum in x 2 space indicates how well these parameters are determined. One 
can also view the dependence of the x 2 m each 7rN charge channel, or the DR pseudo-data, separately to look for 
systematic variations. 

We found the mapping procedure to be necessary for {g 2 /4tt, a^ + ) due to computational instabilities which arise 
occasionally. One could, in principle, fix all the parameters and follow the above procedure over a multi-dimensional 
grid, and tests were done to investigate this, but it was found that the above procedure was sufficient to ensure stable 
fits with well-satisfied DRs in all cases, while keeping the x 2 St to the data optimal. Adding more, or fixing more 
constraints, increases the data fit x 2 ■ 

C. Coulomb Interaction 

The partial-wave and DR analyses are performed using purely isospin-conserving hadronic amplitudes, so Coulomb 
and other isospin-violating contributions must be introduced before constructing observables to compare to the phys- 
ical pion-proton scattering data. There is no unique way to split the interaction into strong and electromagnetic 
parts so in general models must be used. 

In previous analyses, we have used a simple point-source barrier factor correction for all partial waves at all energies, 
as discussed in Ref. j4j. The main isospin-violating contribution is the A ++ — A mass and width difference, which 
we have been treating phenomenologically by splitting the average resonant hadronic P33 partial-wave into and 
i-33 waves for Tt + p and n~p, respectively. The A mass and width differences were found from a best fit to the data. 
A further correction was made to the P33 amplitude to account for the n~p — > channel. The mass (A — A ++ ) 
and width (T° - differences are 1.74±0.15 MeV and 1.09±0.64 MeV, respectively |57| . 

The main criticism of point-source barrier factors is that they exaggerate the S-wave corrections at low energies, 
since the I = partial waves do not vanish at the origin, while the point-Coulomb field diverges there. Consequently, 
in this study, we have used the Nordita [5^] Coulomb prescription, which was used for the Karlsruhe-Helsinki (KH80) 
analyses [5!j. However, the Nordita work supplied corrections for only the largest S- and P-waves up to 530 MeV, 
whereas our analysis extends to 2100 MeV. The Nordita corrections were therefore supplemented with corrections 
calculated from Gibbs extended-source barrier factors [6^ for S- and P- waves at energies above 450 MeV, and D- 
and higher waves at all energies. The Nordita and Gibbs results were smoothly joined in the 450 MeV region. As 
well, an error was found and corrected in the isospin mixing term in one Nordita source |58|. while the pure Coulomb 
part of the (dominant) resonant P33 correction was reduced 15% to take into account intermediate a exchange |61| . 
which was missing in the Nordita work. 

The Nordita corrections already contain effects due to the 7r + — 7r° mass difference. In this scheme, the 7r + channel 
corresponds to the hadronic channel, and all isospin violating effects (due to mass differences and the 771 channel) are 
subsumed in the 7r~ amplitudes. To remain consistent with this scheme, we have redefined our A(1232) resonance 
charge-splitting analogously. This has had the effect of shifting the peak of the hadronic resonance to a slightly lower 
energy. 

Final FA02 Coulomb correction results for the S31 and P33 partial-wave are plotted in Fig. [5] along with the point- 
source barrier factor results. Differences between the FA02 and point-source corrections (not shown) are similarly 
large in the Sn partial-wave, as expected, while for all other waves, the differences are much smaller. However, these 
differences in the SM95 and FA02 correction schemes did not result in large effects on the fits or dispersion relation 
parameters. 
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IV. RESULTS AND DISCUSSION 



A. FA02 versus the SM95 Fit 



The main result of this work has been an energy-dependent solution (FA02) having a \ 2 °f 45874 for 23979 (n^ 
clastic, with pion induced ir° and rj production) data to 2.1 GeV (77N data have been included to 800 MeV). The 
overall x 2 /data was significantly lower than that found in our previously published fit (x 2 /data = 2.2) despite 
the inclusion of additional data. This change is partly a reflection of the database changes mentioned in Section [H] 
Our present and previous energy-dependent solutions are compared in Table [I] As in previous analyses, we have 
used the systematic uncertainty as an overall normalization factor for angular distributions. For total cross sections 
and excitation data, we have combined statistical and systematic uncertainties in quadrature. This renormalization 
freedom allows a significant improvement in our best-fit results, as shown in Table ITT1 In cases where the systematic 
uncertainty varies with angle, this procedure may be considered a first approximation. 

In Table ITTT1 we compare the energy- dependent and single-energy results over the energy bins used in these single- 
energy analyses. Also listed is the number of parameters varied in each single-energy solution. A total of 172 (94 
for isospin 1/2 and 78 for isospin 3/2) parameters were varied in the energy-dependent analysis FA02. The extended 
database has allowed an increase in the number of single-energy values versus our previous result [|| over the same 
energy range. 

Figs. and compare the energy-dependent fits FA02 and SM95 over the full energy region. Partial waves with 
/ < 6 are displayed, whereas the analysis has fitted waves up to I = 6. Deviations from SM95 are visible in amplitudes 
S and P for I = 1/2, and the D waves, mainly near the end point of our analysis. Considering the complicated 
structure associated with an opening ?m channel and N(1535) resonance, the Sn partial-wave shows remarkably little 
deviation from the earlier SM95 fit. 



B. Resonance Parameter Extraction 



The resonance spectrum from our fit has been extracted using three different methods. Poles and zeros have been 
found by continuing into the complex energy plane. These are compiled in Tables IVTl IVIII and IVIIII along with the 
modulus and phase of all pole residues. The complicated interplay of poles and zeros is displayed in Fig. [S] for the 
Sn and Pn partial waves. Clearly, for these partial waves and their associated resonances, a simple Breit-Wigner 
paramcterization should be avoided. An error in the routine used to determine pole residues for SM95 was found |62l | 
and corrected. The most significant change is seen in the A(1232) residue, which now agrees with most previous 
determinations. 

More commonly used, though more model-dependent, are Breit-Wigner parameters for the resonances. Here, as for 
the SM95 fit, a Breit-Wigner plus background contribution was initially fitted to the single-energy amplitudes, over 
varying ranges of energy. However, given the scatter often found in the single-energy values, we have rejected this 
method in favor of another. In Tables IIVI and we give resonance parameters obtained from a Breit-Wigner plus 
background (using a product S-matrix approach, S — SrSb) representation applied directly to the data. Here, Sr = 
1 + 2 i T R , with 

ji _ 

R W R -W -i (T w /2 + T 7 /2) " 1 ' 

The total width is given by T = T n + T/, where 

r w = P7r TR, (5) 



r 7 = Pi r(i-R), 

with R being the branching fraction to 7rN. The background T-matrix is given in terms of a K-matrix, Tg = — 
iKs)^ 1 , with Kb = A + B(W — Wr) + iC. Only one background parameter was necessary for the I = 3/2 waves. 
Data were then fitted using this representation for a particular resonant partial wave. The remaining waves were 
fixed to values found in the full global analysis. Values for the background parameters, energy ranges over which fits 
were performed, and % 2 comparisons are given in Tables llXl and Ixl In Fig. 13 S-, P-, D-, and F-wave resonances are 
compared in terms of their pole positions versus a fitted Breit-Wigner mass and width. 
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In Table llVl results for the N(1535) and N(1650) resonances are significantly different than those reported in the 
SM95 analysis. The widths of both states are now closer to their RPP averages. The N(1650) continues to show a 
nearly elastic behavior in the FA02 solution. This structure also appears to be difficult to parameterize, requiring the 
most background parameters in the present representation. Breit-Wigner parameters are absent for the A(I920) as 
this fit produces a mass about 500 MeV greater than suggested by the pole position. 

In Table Ixll we give resonance parameters determined from a fit to data without DR constraints. In most cases, 
there is reasonable agreement with the constrained fit. An exceptions appears to be the Pi3(1720), which has a 
relatively weak coupling to ttN states, and has an atypical resonance signature. 



C. Charge-Symmetry Violation in tvN Scattering 

The issue of charge symmetry violation (CSV) is fundamental to our understanding of hadronic interactions, and 
many experimental and theoretical studies have addressed this issue (see review [63J). In the framework of QCD, 
CSV arises from the mass difference between the u and d quarks. The other principal cause for CSV comes from 
the electromagnetic interaction, as described in Sec. IIII CI Weinberg [64[ pointed out that the effective chiral 
7rN Lagrangian, coming from QCD, contains a term which violates charge symmetry (see also, a recent review by 
Meissner 65] ). Thus, not only are there kinematic reasons for CSV due to the mass differences within baryon 
multiplets, but direct CSV effects should exist as well. Recent analyzes |6^,|63 of the triangle identity in low-energy 
7rN scattering (below 70 MeV) data have found some indications for significant (6-7%) direct CSV effects in the 
strong-interaction sector [ofil W\. Fettes and Meissner [6^| have found a somewhat smaller violation (~2.5%) in a 
third order chiral perturbation theory calculation without electromagnetic effects, and even smaller violations once 
the latter were included [fjjj. 

Because our formalism does not employ CSV beyond the Coulomb interaction, which we take into account in our 
treatment of the database, we have generated a test fit (X370) of the full database excluding charge-exchange data. 
Differences between FA02 and X370 are therefore signal incompatibilities between the elastic and charge- exchange 
data, which could be an indication of CSV. 

In Table IXHI we compare the FA02 and X370 solutions. For the energy range associated with meson factories 
(below 500 MeV), there is little difference between FA02 and X370. Above 500 MeV, a significant x 2 contribution 
comes from the old Rutherford High Energy Laboratory (RHEL) ir + p P measurements where there are 1976 P data 
from 480 to 2080 MeV [70j favoring by x 2 ~200 the fit X370 versus FA02. Conversely, for charge-exchange da/dtl two 
old RHEL sets (416 cross sections from 900 to 2050 MeV E3. 1 56 backward cross sections from 480 to 870 MeV [37j, 
and one KEK set (72 cross sections from 1830 to 2040 MeV |71|) are more poorly represented, by \ 2 ^1900, in X370 
versus FA02. At very low energies, there exist rather few charge-exchange data, and almost no polarization data, so a 
precise test of the Fettes and Meissner j6^,|6|j prediction is not possible. As a result, we cannot claim any compelling 
evidence for sizable CSV effects in 7rN elastic scattering. This is consistent with recent findings for it2N and ir3N 
systems [72|. 



D. Dispersion Relation Parameters 

Results for the forward C (w), forward derivative E^fa), and fixed-t C (u, t), Hiiper, and B + {v,t) DRs are 
summarized graphically in Figures 121 — 141 ITTJ1 and ^2 respectively. The form of B + (i/,t) DR displayed in Fig. 0] was 
used in Ref. [73 in a determination of the coupling constant. This form was not used as a constraint in FA02, but is 
included here for its utility in illustrating the uncertainty associated with g 2 /4ir. These DRs are well satisfied over 
the whole constraint region up to a kinetic energy T„ = 1 GeV and a four-momentum transfer t = —0.4 (GeV/c) 2 . 
The consistency of FA02 with these constraints is much better than our previous solution SM95 Q and the solutions 
KA84 [54| and KH80 |5fj. Due to the energy-dependent parameterization of our partial waves, all fixed-t DRs are 
well satisfied up to t = —0.5 (GeV/c) 2 , with all isovector DRs reasonably well satisfied up to ~2 GeV i.e. the entire 
data range. The more sensitive isoscalar DRs remain reasonably well satisfied up to 1100—1300 MeV; constraints for 
the forward C+(u>) DR extended to 0.8 GeV. 

Figure E| compares the "Chew-Low" [5(J and scattering length (volume) plots for the FA02 and Karlsruhe [54| S- 
and P-wave amplitudes. Exact agreement with the values given in Figgis not expected, as the method is different. 
The FA02 Chew-Low plots have a very similar behavior to the partial-wave dispersion relation (PWDR) constrained 
KA84 solution. This indicates that FA02 at least approximately satisfies a P-wave PWDR. The present solution fixes 
a long-standing discrepancy in the P13 Chew-Low behavior observed in our previous solutions 0, IM Q ■ 
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The present solution FA02 yields the DR parameters g 2 /47r = 13.75±0.10 for the pion-nucleon coupling constant, 
= 0.133±0.001 n~ 3 (a7 + = -0.074±0.001 ^~ 3 ) for the P-wave scattering volumes, and 3a^ + = 0.2650±0.0014 yT 1 
and a,,-,, = 0.0870±0.0013 [T 1 (implying an, = -0.0010±0.0012 yT 1 ) for the S-wave scattering lengths. The coupling 
constant confirms our earlier result in Ref . [5| and is in line with other determinations 0, FiL The scattering lengths 
are in agreement with those of the PSI pionic- hydrogen and pionic-deuterium experiments |55| . 

Our DR procedure constrains the 7r _ p scattering length to agree (within some error bound) with the result from 
the PSI pionic-atom experiment 55] (a^-p = 0.0883±0.0008 /i^ 1 ), so the above agreement in a T ~ p is not surprising. 
However, a solution constructed to be otherwise identical to FA02, but where not constrained to the 

PSI value, yielded a best fit with a^-p = 0.0856±0.0010 fi , indicating the value sought by the scattering data 
alone. Consequently, a compromise value of ~0.087 a little lower than the PSI result was chosen as the constraint 
value for FA02. Ericson, Loiseau, and Wycech have published [Tfj a reanalysis of the PSI result [5^], and have 
obtained a^ = 0.0883±0.0027 /i -1 and a^- p = 0.0870±0.0005 /i -1 , the latter reasonably consistent with our result 
0.0856±0.0010 fi- 1 using data alone. 

The systematic check described above is one of many that we have applied to our solution. A large number of test 
solutions were obtained with varying fitting and constraint conditions, to investigate systematic uncertainties in the 
extracted DR parameters. These systematics checks are discussed in detail below. 

• A solution (P370) using the point-Coulomb corrections employed in our previous solutions Refs. [H, H> El gave 
an overall \ 2 fit very similar to FA02. Differences in the resulting parameters were not large. For g 2 /Air and 
a^, , this solution gave 13.69 and 0.0867 M -1 , respectively. Despite seemingly large differences between the 
point-Coulomb and Nordita+Gibbs jH, |6(| schemes (see Fig. |SJ), the consistency of results indicates that our 
Coulomb correction scheme is adequate to perform this type of partial- wave analysis. 

• In an isospin-invariant framework, charge-exchange data are not required to reconstruct the amplitudes. As 
discussed in Sec. IIV CI a solution was constructed (X370) after charge-exchange data were removed from the 
fitted database. Despite this drastic change, the resulting DR parameters were quite consistent with those of 
FA02. 

• A less drastic change in the database was investigated in the region of the A(1232) resonance. The issue of 
discrepant data sets on the rising part of the resonance was discussed in Sec.^O To quantify this effect, solutions 
were constructed employing the data sets of Pedroni et al. 29] , Brack et al. Il4l |3(j , and Pavan et al. , while 
floating (applying a normalization error of 100%) those of Bussey et al. |32j ] and Carter et al. |3l|. and vice 
versa. The 7r~p scattering-length constraint was held constant in both cases. The resulting solutions differed in 
the coupling constant by Ag 2 /47r = 0.07, while the scattering lengths remained consistent. Since FA02 fits all 
the aforementioned data sets together, an uncertainty of ±0.04 was ascribed to this systematic effect. 

• The DR constraints for the FA02 solution were applied up to 1 GeV, and from t = to —0.4 (GeV/c) 2 . Solutions 
constructed with a reduced range of 25 to 600 MeV and t to —0.3 (GeV/c) 2 again yielded results consistent 
with FA02, the latter ranges having been employed in solution SM95 [5j. 

• Dispersion relation constraints were both strengthened and weakened by a factor of 2, yielding results consistent 
with FA02. It should be noted, however, that constraining too tightly yields a very large data \ 2 an d numerically 
unstable fits. 

• In addition to the above test, a solution was constructed employing no DR constraints. Here, the fit to data was 
naturally much better than in FA02, but the improvement was less dramatic (A\ 2 = 1036 for 21808 data points) 
than one might expect, with little difference obvious in a plot of the partial waves. The numerically sensitive C + 
and E + isoscalar dispersion relations were understandably not well satisfied; however, the less sensitive C~ and 
E~ isovector dispersion relations remained rather well satisfied over the full energy range, as were the Hiiper 
and B + (j/,t) relations up to ~800 MeV and t = —0.3 (GeV/c) 2 . The extracted constants were ag + = 0.086 /x _1 , 
aj~ + = -0.073 fi- 3 , and g 2 /47r = 13.62. 

• The D-waves below ^250 MeV are too small to be accurately determined from fits to data alone. The forward 
derivative E* DRs have an enhanced sensitivity to higher partial waves (~ ^ 3 ) and so help to constrain the D- 
waves. To test the dependence of the DR parameters on these low energy D-waves, a solution was constructed 
identical to FA02 except constrained rigorously (1%, compared to 5% for FA02) to the low energy Koch D-wave 
results 0. The DR parameters changed minimally, the largest changes being 0.002 for both Aa^ + , and Aa° + . 
However, the data \ 2 was increased both above 400 MeV (A\ 2 = ~600) and below by ^200, with an increase 
of ^300 in the fit to DR pseudo-data constraints. The fit thus exhibits a clear preference for the FA02 D-waves. 



9 



The central values and uncertainties of the DR parameters have been estimated, taking into consideration the 
systematic checks outlined above, and from other checks within the FA02 solution itself. These include differences in 
parameter values at the \ 2 minima for each of the charge channels and the DR pseudo-data. For the coupling constant 
g 2 /4ir, we also considered differences in the value extracted from the Hiiper (Tig. llU|) and B + (y, t) fFig. lllfl dispersion 
relations, which, along with the database changes around the A(1232), were the dominant sources of uncertainty. 
The statistical uncertainty stemming from A\ 2 = 1 variations from the minima is negliglible compared to systematic 
effects. 

These tests for systematic effects improve our confidence in the values of parameters extracted from FA02. It should 
be noted, however, that the scattering lengths are tied to the PSI pionic-atom value for the 7r~p scattering length. 
A new more precise experiment is planned and if the result were to change several standard deviations from 
the present value, this would necessitate a re-analysis. Nonetheless, the value of the pion-nucleon coupling constant 
would not be expected to change, as it is observed to be robust with respect to modest changes in the 7r~p scattering 
length. 

E. The n~p — > rjn Channel 

Our fit to a representative set of ir~p — > r\n cross sections is displayed in Fig. 1131 In general, over this region, 
our results are qualitatively similar to other recent multi-channel analyses which have included additional reaction 
data [3 • Compared to the 7rN elastic result, the overall x 2 f° r this channel is slightly higher. However, we feel this is 
largely due to problems in the database. 

Our coupled-channel a ppr oach allows the determination of a number of amplitudes related to the 77N interaction. 
Fig. 1141 gives an Argand |78j plot representation of the energy-dependent fit FA02 over the energy region including 
tt~p — ► -qn data. Fig. IT31 gives a more detailed view of our 7rN elastic and irN — ► r/N results for the Su and Di 3 
partial waves. 

V. SUMMARY AND CONCLUSION 

We have fitted the existing 7rN elastic scattering database, employing a complete set of DR constraints, up to = 
1 GeV and t = —0.4 (GeV/c) 2 . Data from the reaction ir~p — > rjn have been included for the first time in an analysis 
of this type. These improvements have allowed us to more carefully examine the N(1535), which is nearly obscured (in 
the 7rN elastic reaction) by the opening of the r]N threshold. Remarkably, the resulting Su partial wave shows little 
change over this energy region. However, the N(1535) resonance width has changed (increased) dramatically, this 
being due mainly to our new method for fitting Breit-Wigner parameters. This result should be taken into account 
by any multi-channel analysis which fits single-energy partial-wave amplitudes. 

Our fits without DR constraints and without charge- exchange data have also yielded interesting results. In most 
cases, an extensive use of DRs had little effect on the extracted resonance spectrum. For weak or non-canonical 
structures, however, dispersion relations may play a more important role. 

In the absence of charge symmetry violation (CSV), beyond Coulomb effects, one should be able to predict the 
charge-exchange observables based on a fit to elastic scattering data. In the low-energy region, where CSV effects are 
expected to be most important, our fit excluding charge-exchange data is reasonably predictive, and could be useful 
in a more refined study of this issue. 

An extensive list of tests designed to check for systematic effects (changes to the database, Coulomb-correction 
scheme) has (a) revealed the DR parameters to be fairly robust with respect to these effects, and (b) suggested a way 
to quantify these systematic uncertainties. As the FA02 the solution exhibits good consistency with the complete set 
of forward and fixed-t DRs, we have added confidence in our results. In particular, the 7rNN coupling has remained 
quite stable, changing little from our SM95 determination. Our DR analysis and systematic checks will next be 
applied to an extraction of the nucleon cr-term (see, e.g. Ref. H^|). The extraction of this quantity is sensitive to the 
fine details of the analysis, and a careful examination of the cr-term is in preparation. 
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TABLE I: Comparison of present (FA02) and previous (SM95 0, FA93 H, SM90 @], and FA84 Q) energy-dependent partial- 
wave analyses of elastic Ap, charge-exchange (7r°n), and 7r~p — > 7771 (rjn) scattering data. For FA02 solution, r/N data has 
been included to 800 MeV. The \ 2 values for the previous solutions correspond to our published results. N prm is the number 
of parameters (7 = 1/2 and 3/2) varied in the fit. 

Solution Range (MeV) X 2 / 7r+ P X 2 / 7r_ P X 2 / 71 " ™ X 2 /v n N prm 



FA02 


2100 


21735/10468 18932/9650 4136/1690 


439/173 94/78 


SM95 


2100 


23593/10197 18855/9421 4442/1625 


94/80 


FA93 


2100 


23552/10106 20747/9304 4834/1668 


83/77 


SM90 


2100 


24897/10031 24293/9344 10814/2132 


76/68 


FA84 


1100 


7416/ 3771 10658/4942 2062/ 717 


64/57 



TABLE II: Comparison of x 2 /data for normalized (Norm) and unnormalized (Unnorm) data used in the FA02 solution. 



Reaction Norm Unnorm 



TT + p — > 


n + p 


2.1 


9.3 


7T~p — > 


n~p 


2.0 


7.1 


7T~p — ► 


A 


2.4 


9.5 






2.5 


4.6 
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TABLE III: Single-energy (binned) fits of combined elastic tt p, charge- exchange, n~p — > r/n scattering data, and \ 2 values. 
Np r m is the number of parameters varied in the single-energy fits, and \\ is given by the energy-dependent fit, FA02, over the 
same energy interval. 



TV (MeV) Rang. 


i (MeV) N prm 


X 2 /data Xe 


TV (MeV) Range (MeV) N prm 


X 2 /data 


xl 


30 


26 


- 33 


4 


171/141 


218 


795 


793 


- 796 


30 


204/165 


301 


47 


45 


- 49 


4 


74/ 82 


80 


820 


813 


- 827 


30 


399/302 


436 


66 


61 


- 69 


4 


183/136 


198 


868 


864 


- 870 


31 


294/211 


389 


90 


87 


- 92 


4 


116/111 


154 


888 


886 


- 890 


32 


174/144 


291 


112 


107 


- 117 


6 


114/ 93 


115 


902 


899 


- 905 


34 


539/416 


891 


124 


121 


- 126 


6 


84/ 60 


92 


927 


923 


- 930 


34 


234/200 


378 


142 


139 


- 146 


6 


230/159 


213 


962 


953 


- 971 


34 


387/299 


593 


170 


165 


- 174 


6 


177/141 


170 


1000 


989 


-1015 


36 


691/423 


847 


193 


191 


- 195 


6 


103/107 


119 


1030 


1022 


-1039 


38 


286/272 


396 


217 


214 


- 220 


6 


116/109 


143 


1044 


1039 


-1048 


38 


362/243 


520 


238 


235 


- 241 


6 


124/115 


149 


1076 


1073 


-1078 


38 


240/218 


432 


266 


263 


- 271 


6 


162/123 


175 


1102 


1099 


-1103 


39 


232/173 


346 


292 


291 


- 293 


8 


157/129 


199 


1149 


1147 


-1150 


40 


339/210 


468 


309 


306 


- 310 


8 


160/140 


174 


1178 


1165 


-1192 


41 


759/394 


925 


334 


332 


- 335 


8 


94/ 58 


130 


1210 


1203 


-1216 


42 


286/233 


364 


352 


351 


- 352 


9 


82/110 


99 


1243 


1237 


-1248 


44 


473/283 


586 


389 


387 


- 390 


10 


31/ 28 


74 


1321 


1304 


-1337 


44 


731/401 


946 


425 


424 


- 425 


10 


151/139 


191 


1373 


1371 


-1374 


46 


331/166 


613 


465 


462 


- 467 


13 


358/120 


451 


1403 


1389 


-1417 


48 


544/408 


811 


500 


499 


- 501 


16 


161/136 


209 


1458 


1455 


-1460 


48 


279/258 


504 


518 


515 


- 520 


18 


104/ 78 


144 


1476 


1466 


-1486 


48 


510/344 


682 


534 


531 


- 535 


18 


131/128 


184 


1570 


1554 


-1586 


48 


834/546 


1062 


560 


557 


- 561 


19 


310/151 


582 


1591 


1575 


-1606 


48 


402/336 


588 


580 


572 


- 590 


20 


382/286 


556 


1660 


1645 


-1673 


48 


510/391 


732 


599 


597 


- 600 


22 


253/151 


382 


1720 


1705 


-1734 


48 


388/279 


510 


625 


622 


- 628 


24 


125/ 94 


204 


1753 


1739 


-1766 


48 


654/439 


864 


662 


648 


- 675 


25 


554/352 


735 


1838 


1829 


-1845 


48 


444/290 


698 


721 


717 


- 725 


28 


160/129 


217 


1875 


1852 


-1897 


48 


948/674 1284 


745 


743 


- 746 


28 


159/100 


286 


1929 


1914 


-1942 


48 


857/501 1187 


765 


762 


- 767 


30 


190/167 


289 


1970 


1962 


-1978 


48 


461/271 


684 


776 


774 


- 778 


30 


220/155 


302 


2026 


2014 


-2037 


48 


368/320 


652 
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TABLE IV: Resonance couplings from a Breit-Wigner fit to the FA02 solution [GW], our previous solution SM95 [VPI] and 
an average from the Review of Particle Properties [RPP] (in square brackets). Masses Wh, half-widths r/2, and partial 
widths for r,rjv/r are listed for isospin 1/2 baryon resonances A*. 



R.csoii3,nc6 


Wb 


r/2 


r^w/r 

7TiV / L 


Ref 




(MeV) 


(MeV) 






P„ flAAfW 


140o.Uzt4.0 


iou_tio 


n 7p;n-i-0 09/1 


n w 
v i vv 




1467 


220 


0.68 


VPI 




1 J,*fc«±U| 




To fi^l 


RPP 
xv l i 


n, ^ (\ ^90^ 


lOlD.oztU.o 


'xy.oznl.o 


n A/io-i-0 oo^ 


P W 

Vjr VV 








fil 

U.U1 


VPT 

V 1 1 




[10ZUJ 


ffinl 

[OUJ 


To ncil 


RPP 


a, . (\ 

^lOOO ) 


1 Pl/i « 7-1- O 9 


oy .UznO.o 


n Q^ri-l-0 OOQ 


P W 

Vjr VV 




1UUU 


~X\ 

UJ 


^1 


VPT 

V 1 1 




n ^^i 

[luouj 




TO A K\ 


RPP 

XV x^ 


Oil I, 1D0U ) 


1 RPil 9-1-/1 7 
lD01.Zzt4. / 


OD.o_to.o 


1 000 

1 .uuu 


P W 
Vjr VV 




1667 


45 




VPI 




[1UUUJ 


[751 


\C\ 791 
[U. (Z\ 


RPP 


n, ^ f i ^7^ 


ID ( D.ZztU.D 


7c q_i_i c: 


n /loo-i-o 009 


p w 

Vjr VV 




1673 


77 


0.38 


VPI 




[1U I UJ 


['°J 


[0 4^1 


RPP 


"15 ^lOoUj 


i coo 9-|_n 7 


£7 9-1- 1 Q 


n ^70-1-0 oo/i 

U.u f UztU.UU^ 


P W 
vjr VV 




1UI o 


uo 


u.uo 


VPT 

V i X 




[1680] 


[65] 


[0.65] 


RPP 


Pia(1720) 


1749.6±4.5 


128±11 


0.190±0.004 


GW 




1820 


177 


0.16 


VPI 




[1720] 


[75] 


[0.15] 


RPP 


G ir (2190) 


2192.1±8.7 


363±31 


0.230±0.002 


GW 




2131 


238 


0.23 


VPI 




[2190] 


[225] 


[0.15] 


RPP 


ffi 9 (2220) 


2270±11 


183±21 


0.200±0.006 


GW 




2258 


167 


0.26 


VPI 




[2220] 


[200] 


[0.15] 


RPP 


Gi 9 (2250) 


2376±43 


462±89 


0.110±0.004 


GW 




2291 


386 


0.10 


VPI 




[2250] 


[200] 


[0.10] 


RPP 



TABLE V: Parameters for isospin 3/2 baryon resonances A*. Notation as in Tabic II VI 



Resonance 


W R 


r/2 


r^jv/r 


Ref 




(MeV) 


(MeV) 






P 33 (1232) 


1232. 9±1. 2 


59. Oil. 1 


1.000 


GW 




1233 


57 


«1.0 


VPI 




[1232] 


[60] 


[>0.99] 


RPP 


531(1620) 


1614.1±1.1 


70.5±3.0 


0.310±0.004 


GW 




1617 


54 


0.29 


VPI 




[1620] 


[75] 


[0.25] 


RPP 


D 33 (1700) 


1687.9±2.5 


182.4±8.3 


0.150±0.001 


GW 




1680 


136 


0.16 


VPI 




[1700] 


[150] 


[0.15] 


RPP 


F 35 (1905) 


1855.7±4.2 


167±11 


0.120±0.002 


GW 




1850 


147 


0.12 


VPI 




[1905] 


[175] 


[0.10] 


RPP 


D 35 (1930) 


2046±45 


201±99 


0.040±0.014 


GW 




2056 


295 


0.11 


VPI 




[1930] 


[175] 


[0.15] 


RPP 


F 37 (1950) 


1923.3±0.5 


139.1±1.5 


0.480±0.002 


GW 




1921 


116 


0.49 


VPI 




[1950] 


[150] 


[0.37] 


RPP 
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TABLE VI: Pole positions from the solution FA02 [GW], our previous solution SM95 [VPI] @, and an average from the Particle 
Data Group [RPP] Q (in square brackets). Re (Wr) and Im (Wi) parts are listed for isospin 1/2 baryon resonances. Second 
sheet poles are labeled by f. Modulus and phase values are listed for the nN elastic pole residue. 



Wave 


Wr 


-Wr 


Modulus Phase 


Ref 


Wave 


Wr 


-Wi 


Modulus Phase 


Ref 




(MeV) (MeV) 


(MeV) 


(dee) 






(MeV) (MeV) 


(MeV) 


(dee) 




5xi 


1526 


65 


33 


14 


GW 




1659 


73 


29 


-22 


GW 




1501 


62 


31 


-12 


VPI 




1663 


76 


29 


-6 


VPI 




[1505] 


[85] 

L J 






RPP 




[1660] 


[701 






RPP 


5ii 


1653 


91 


69 


-55 


GW 


Fl5 


1678 


60 


43 


1 


GW 




1673 


41 


22 


29 


VPI 




1670 


60 


40 


1 


VPI 




[1660] 


[80] 






RPP 




[1670] 


[60] 






RPP 


Pa 


1357 


80 


36 


-102 


GW 


Fit 


1779 


124 


47 


-61 


GW 




1346 


88 


42 


-101 


VPI 




1793 


94 


27 


-56 


VPI 




[1365] 


[105] 






RPP 












RPP 


P" 


1385 


83 


82 


-51 


GW 


G17 


2076 


251 


68 


-32 


GW 




1383 


105 


92 


-54 


VPI 




2030 


230 


46 


-23 


VPI 












RPP 




[2050] 


[225] 






RPP 


Pl3 


1655 


139 


20 


-88 


GW 


Gig 


2238 


268 


33 


-25 


GW 




1717 


194 


39 


-70 


VPI 




2087 


340 


24 


-44 


VPI 




[1700] 


[125] 






RPP 




[2140] 


[240] 






RPP 


Dl S 


1514 


51 


35 


-6 


GW 


Hig 


2209 


282 


96 


-71 


GW 




1515 


55 


34 


7 


VPI 




2203 


268 


68 


-43 


VPI 




[1510] 


[57] 






RPP 




[2170] 


[235] 






RPP 



TABLE VII: Pole positions for isospin 3/2 baryon resonances A*. Notation as in Table IVT1 



Wave 


W R 


-Wi 


Modulus Phase 


Ref 




(MeV) (MeV) 


(MeV) 


(deg) 




031 


1594 


59 


17 


-104 


GW 




1585 


52 


14 


-121 


VPI 




[1600] 


[57] 






RPP 


P31 


1748 


262 


48 


158 


GW 




1810 


247 


53 


-176 


VPI 




[1855] 


[175] 


— 




RPP 


P33 


1210 


50 


53 


-47 


GW 




1211 


50 


38 


-22 


VPI 




[1210] 


[50] 






RPP 


D33 


1617 


113 


16 


-47 


GW 




1655 


121 


16 


-12 


VPI 




[1660] 


[100] 






RPP 


D35 


1966 


182 


16 


-21 


GW 




1913 


123 


8 


-47 


VPI 




[1890] 


[125] 






RPP 


F35 


1825 


135 


16 


-25 


GW 




1832 


127 


12 


-4 


VPI 




[1830] 


[140] 






RPP 


F37 


1874 


118 


57 


-34 


GW 




1880 


118 


54 


-17 


VPI 




[1885] 


[120] 






RPP 
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TABLE VIII: Zero positions from the solution FA02 [GW] and our previous solution SM95 [VPI] [Ej. Re (Wr) and Im (Wi) 
parts are listed for isospin 1/2 and 3/2 baryon resonances. 



VV cLVc 


Wr 


-Wi 


Ref 




(MeV) (MeV) 




OH 


1578 


38 


V 1 vv 




1582 


54 


VPT 
v jt 1 




1695 


43 


VPI 


Pl3 


1585 


51 


GW 




1618 


85 


VPI 


D 13 


1759 


64 


GW 




1751 


87 


VPI 


Fib 


1765 


66 


GW 




1775 


57 


VPI 


S31 


1580 


36 


GW 




1579 


30 


VPI 


P31 


1826 


197 


GW 




1863 


170 


VPI 


D 35 


1806 


107 


GW 




1827 


69 


VPI 



21 



TABLE IX: Comparison of FA02 and Breit-Wigner+background fits. Background parameters for isospin 1/2 baryon resonance 
fits (see text and associated Table If VI . "Data" refers to the number of scattering data used in the fit. 



Resonance 


Wmin 


Wmax BW fit 


FA02 


Data 


A 


B 


C 




(MeV) (MeV) 


x 2 


x 2 










Pii(1440) 


1350 


1550 


5556 


5587 


2393 


-0.270±0.026 






Di 3 (1520) 


1480 


1560 


3238 


3341 


1448 


-0.035±0.013 






Su(1535) 


1490 


1590 


3566 


3657 


1632 


0.342±0.018 


4.530±0.851 




5n(1650) 


1600 


1720 


6992 


7218 


3394 


-0.584±0.073 


2.415±1.356 


0.414±0.024 


D 1S (1675) 


1610 


1730 


7343 


7359 


3546 


-0.034±0.004 






Fi 5 (1680) 


1620 


1730 


6683 


6650 


3172 


O.OlliO.011 


1.192±0.575 




Pis(1720) 


1650 


1790 


8321 


8484 


4160 


-0.127±0.008 


1.463±0.239 




Gi 7 (2190) 


2049 


2249 


8017 


8040 


3774 


0.024±0.009 






#19(2220) 


2140 


2250 


5675 


5682 


2554 


0.028±0.010 






Gi 9 (2250) 


2010 


2258 


11020 


11025 


5330 


-0.046±0.010 







TABLE X: Comparison of FA02 and Breit-Wigner+background fits. Background parameters for isospin 3/2 baryon resonance 
fits (see text and associated Table 0. "Data" refers to the number of scattering data used in the fit. 



Resonance 


Wmin 


Wmax BW fit 


FA02 


Data 


A 




(MeV) (MeV) 


x 2 


x 2 






P 33 (1232) 


1180 


1270 


1180 


1185 


920 


0.035±0.017 


S* 3 i(1620) 


1570 


1680 


5187 


5212 


2321 


-0.851±0.013 


£>33(1700) 


1610 


1770 


9624 


9690 


4725 


-0.145±0.003 


F 35 (1905) 


1770 


1920 


8069 


8096 


3791 


-0.092±0.004 


-D 35 (1930) 


1870 


2100 


9912 


9881 


5059 


-0.063±0.014 


F 37 (1950) 


1800 


2000 


10623 


10552 


4951 


0.024±0.006 
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TABLE XI: Parameters for low-lying resonances of isospin 1/2 and 3/2 in a fit unconstrained by dispersion relations. 

Resonance Wr F/2 F^n/T 
(MeV) (MeV) 



Pn(1440) 1473.0 200.9 0.695 

Di 3 (1520) 1516.2 50.8 0.655 

5n(1535) 1545.4 87.0 0.385 

5n(1650) 1658.9 55.3 1.000 

Di B (1675) 1673.7 75.9 0.396 

Fi 5 (1680) 1677.8 64.1 0.692 

Pi 3 (1720) 1807.1 239 0.195 

P 33 (1232) 1233.3 59.3 1.000 

5 3 i(1620) 1614.0 70.1 0.316 

D 33 (1700) 1684.9 168.6 0.151 

F 35 (1905) 1861.3 159 0.119 



TABLE XII: Comparison of x 2 for the FA02 (full data set) and X370 (no charge-exchange data) solutions to 2.1 GeV (iv~p — > r\n 
data to 800 MeV). Results to 500 MeV have been listed in brackets, where one observes little or no difference in the fit quality 
between FA02 ands X370, even for the charge-exchange database. "Data" refers to the number of scattering data used in the 
fit. 



Reaction 


Observable 


FA02 


X370 


Data 


ir + p — > ir + p 


da/dfl 


14574 (2043) 14624 (2033) 7246 (977) 




O 'tot 


257 (115) 


220 (115) 


105 (60) 




P 


6801 (533) 


6629 (536) 


3012 (495) 




R 


28 (20) 


29 (20) 


48 (26) 




A 


39 (24) 


42 (25) 


48 (26) 


-K~p — ► ir~p 


da/dQ 


13799 (1443) 13760 (1453) 7331 (773) 




Otot 


531 (111) 


445 (110) 


151 (59) 




P 


4441 (562) 


4501 (564) 


2045 (337) 




R 


107 (24) 


104 (23) 


61 (27) 




A 


49 (14) 


47 (14) 


60 (26) 


iv~p — > 7T°n 


dti/dQ 


3180 (697) 


5093 (739) 


1333 (379) 




a tot 


24 (23) 


24 (23) 


34 (33) 




P 


933 (194) 


1175 (215) 


323 (159) 


ir~p — » ijn 


da/dQ 


373 


376 


138 




a tot 


67 


75 


35 
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FIG. 1: Energy-angle distribution of recent (post-1995) data: (a) -K + p and (b) 7v p clastic scattering, (c) charge-exchange 
ir~p — > 7r°n, and iv~p — > rjn (all available). Total cross sections are plotted at zero degrees. 
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FIG. 2: Subtraction constants as a function of pion kinetic energy for the forward C dispersion relations plotted in terms of 
the S-wave scattering lengths , and other combinations, in inverse pion mass units. Horizontal lines represent the least-square 
averages of individual values. 
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FIG. 3: Subtraction constants as a function of pion kinetic energy at four values of four-momentum transfer t [in (GeV/c) 2 } for 
the fixed-t (a) C + (u, t) and (b) C~(v, t) dispersion relations. Horizontal lines represent the least-square averages of individual 
values. 
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FIG. 4: Forward derivative (a) E + (ui) and (b) E~(ui) dispersion relations. ReE (ImE) are plotted by solid (dashed) lines, the 
principal value integral, PVi, by dash-dotted lines, and the pole term by short-dash-dotted lines. The respective subtraction 
constants are shown as horizontal solid lines. 
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FIG. 5: Comparison of the point-Coulomb (SM95) and Nordita-Gibbs (FA02) charge corrections A = 5 — 5h, the subscript H 
denoting an hadronic phase, for the (a) 5*3i and (b) P33 partial waves. In both plots, dash-dotted and short-dash-dotted (solid 
and dotted) curves give point-Coulomb 5] (FA02) corrections for ir + p and iv~p, respectively. 
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FIG. 6: Isospin 1/2 partial-wave amplitudes {hn,2j) from T w = to 2.1 GeV. Solid (dashed) curves give the real (imaginary) 
parts of amplitudes corresponding to the FA02 solution. The real (imaginary) parts of single-energy solutions are plotted 
as filled (open) circles. The previous SM95 solution [f| is plotted with long dash-dotted (real part) and short dash-dotted 
(imaginary part) lines. The dotted curve gives the unitarity limit (ImT — T*T) from FA02. All amplitudes are dimensionless. 
Vertical arrows indicate Wr and horizontal bars show full T/2 and partial widths for T^n associated with the FA02 results 
presented in Table |TV| (a) Su, (b) P lly (c) P 13 , (d) D 13 , (e) D 15 , (f) Fig, (g) Gvr, and (h) G w . 
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FIG. 7: Isospin 3/2 partial-wave amplitudes (F21.2.7) from = to 2.1 GeV. Notation as in Fig. © Vertical arrows indicate 
Wr and horizontal bars show full F/2 and partial widths for F^jv associated with the FA02 results presented in Fable [V] (a) 
Sax, (b) P31, (c) P33, (d) D33, (e) D35, (f) F35, and (g) F 37 - 
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FIG. 8: Contour plot of ln\T\ displaying complex plane structures, (a) Sn. Poles: Wp — 1526 - i65 MeV and Wp = 1653 - 
i91 MeV, shown as stars. Zero: Wz = 1578 - i38 MeV, represented as an open triangle; (b) Pn, first sheet. Pole: Wp = 1357 
- i80 MeV, shown as a star; and (c) Pn, second sheet. Pole: Wp = 1385 - i83 MeV, shown as a star. The branch point for 
Pn, 1349 - i50 MeV, is represented as a solid triangle. The 7rA branching cut in (b) and (c) is shown as a solid line. 
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FIG. 9: Comparison of complex plane (bottom panel) and Breit-Wigner (top panel) fits for resonances found in the FA02 
solution. Plotted are the result for (a) S- and P-wave resonances and (b) D- and F-wave resonances. Complex plane poles are 
shown as stars (the boxed star denotes a second-sheet pole). Wr and Wi give real and imaginary parts of the center-of-mass 
energy. The full (nN partial) widths are denoted by narrow (wide) bars for each resonance. 
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FIG. 10: The Hiiper dispersion relation [y] plotted for several values of four-momentum transfer from to —0.4 (GeV/c) 
The y-intercept yields the coupling constant g 2 /Ait. Lines represent the least-square averages of individual values. 
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FIG. 11: Deviations from the mean value of g 2 /An for FA02 from the B + (u, t) DR evaluated over a grid of and 4-momentum 
transfer values. Dash-dotted lines show kinematical limits. 
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FIG. 12: Plots of (a) q 2l+1 cot 8/{s - M 2 ) and (b) ReT/q 2l+1 for FA02 (solid) and KA84 54] (dot-dashed) S- and P-wave 
amplitudes. Scattering lengths (volume) in right panel given in fi~ x ()J,~ 3 ) units. Curves are fits quadratic in TV 




0.3 



70 110 
8 (deg) 



150 



£0.2 



a 

^o.i 

b 

0.0 



601 MeV 




(b) 



30 



70 110 

e (deg) 



150 



0.3 



0.3 




70 110 

6 (deg) 



70 110 

d (deg) 



150 



FIG. 13: Differential cross sections for ir~p -> r/n. (a) 575 MeV, (b) 601 MeV, (c) 612 MeV, and (d) 620 MeV. Experimental 
data are from |34l measurements (filled circles), |3ql (open circles), and |3S|| (open triangles). Solid line shows the FA02 
results. 





FIG. 15: Low-lying states related to the r/N interaction. Plotted are the results for (a) Su and (b) D13 amplitudes. Solid 
(dashed) curves give the real (imaginary) parts of nN — > 7riV amplitudes, dash-dotted (dotted) curves represent the real 
(imaginary) parts of nN — > r/N amplitudes corresponding to the FA02 solution. All amplitudes are dimensionless as in Fig. |S| 
Vertical arrows indicate Wr and horizontal bars show full F/2 and partial widths for T^n associated with the FA02 results. 



